Partonic effects on higher-order anisotropic flows in relativistic heavy-ion collisions 
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Higher-order anisotropic flows 774 and va of charged hadrons in heavy-ion collisions at the Rela- 
tivistic Heavy Ion Collider are studied in a multiphase transport model that has previously been 
used successfully for describing the elliptic flow V2 of identified hadrons in these collisions. We find 
that the same parton scattering cross section of about 10 mb used in explaining the measured V2 of 
charged hadrons can also reproduce the recent data on their V4 and ve from Au + Au collisions at 
y/s — 200 AGeV. It is further found that V4 is a more sensitive probe of the initial partonic dynamics 
in these collisions than vi- Moreover, higher-order parton anisotropic flows are non-negligible and 



satisfy the scaling relation v n ,q(pT) 



n/2 



(Pt), which leads naturally to the observed similar scal- 



ing relation among hadron anisotropic flows when the coalescence model is used to describe hadron 
production from the partonic matter. 

PACS numbers: 25.75.Ld, 24.10.Lx 
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I. INTRODUCTION 

There have been extensive studies on the azimuthal 
anisotropy of hadron momentum distributions in the 
transverse plane perpendicular to the beam direction, 
particularly the lower-order elliptic flow V2, in heavy- 
ion collisions at various energies pj. The hadron trans- 
verse momentum anisotropy is generated from the pres- 
sure anisotropy in the initial compressed matter formed 
in non-central heavy ion collisions [2, and is thus sensi- 
tive to the properties of produced matter. This sensitiv- 
ity not only exists in the larger elliptic flow 0,IEIH0>ISE1 
but also in the smaller higher-order anisotropic flows 
|ToL Ull. lil^. Recently, higher-order anisotropic flows D4, 
Vq, and vs in heavy ion collisions at the Relativistic Heavy 
Ion Collider (RHIC) were studied in more detail using the 
hydrodynamical model |l3[ . and it was suggested that 
i>4 could also provide important information on the col- 
lision dynamics during the early stage. Although the 
hadron anisotropic flows become smaller with increasing 
orders, recent experimental results from STAR Collabo- 
ration [bH have demonstrated that they are still measur- 
able. 

In the present work, the anisotropic flows V2, «4, and 
vq of both partons and hadrons in heavy-ion collisions 
at RHIC are studied in a multiphase transport (AMPT) 
model that includes initial partonic and final hadronic 
interactions Il5l Using the version with string melt- 
ing, i.e., allowing hadrons that are expected to be formed 
from initial strings to convert to their valence quarks and 
antiquarks 0, 0, 113 1 we nn d that the measured data 
from Au + Au collisions at y/s — 200 AGeV can be ap- 
proximately reproduced if a parton scattering cross sec- 
tion of about 10 mb is used in the model. Our results 



show that the hadron U4 is indeed a more sensitive probe 
of the initial partonic dynamics in relativistic heavy-ion 
collisions than 1)2- We also find that higher-order par- 
ton anisotropic flows are non-negligible and satisfy the 
scaling relation v rhq (pT) ~ v 2,qipT), which then leads 
to a similar scaling relation among hadron anisotropic 
flows when hadron production from the partonic matter 
is described by the quark coalescence model. 



II. THE AMPT MODEL WITH STRING 
MELTINGS 

The AMPT model [H |H E3 is a hybrid model 
that uses minijet partons from hard processes and strings 
from soft processes in the heavy ion jet interaction gen- 
erator (HIJING) model j^f as the initial conditions for 
modeling heavy-ion collisions at ultra-relativistic ener- 
gies. Since the initial energy density in Au + Au colli- 
sions at RHIC is much larger than the critical energy den- 
sity at which the hadronic matter to quark-gluon plasma 
transition would occur [2(], Hj| , we shall use the version 
which allows the melting of initial excited strings into 
partons |17| . In this version, hadrons, which would have 
been produced from the HIJING model, are converted 
instead to valence quarks and/or antiquarks. Interac- 
tions among these partons are described by the Zhang's 
parton cascade (ZPC) parton cascade model [2^]. At 
present, this model includes only parton-parton elastic 
scatterings with the in-medium cross section 
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where a 8 = 0.47 is the strong coupling constant, and s 
and t are the usual Mandelstam variables. The effective 
screening mass \x depends on the temperature and den- 
sity of the partonic matter but is taken as a parameter 
in ZPC for fixing the magnitude and angular distribu- 
tion of parton scattering cross section. Since there are 
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no inelastic scatterings, only quarks and antiquarks from 
the melted strings are present in the partonic matter. 
As described in Ref . 0] , the transition from the partonic 
matter to the hadronic matter is achieved using a sim- 
ple coalescence model, which combines two nearest quark 
and antiquark into mesons and three nearest quarks or 
antiquarks into baryons or anti-baryons that are close to 
the invariant mass of these partons. The present coa- 
lescence model is thus somewhat different from the ones 
recently studied extensively [2(| HU , which are used 
only for producing hadrons with intermediate transverse 
momenta. The final-state hadronic scatterings are then 
modeled by a relativistic transport (ART) model |29| . 
Using parton scattering cross sections of 6-10 mb, the 
extended AMPT model is able to reproduce both the cen- 
trality and transverse momentum (below 2 GeV/c) de- 
pendence of the elliptic flow 17] and pion interferometry 
[18f measured in Au+Au collisions at y/s = 130 AGeV 
at RHIC 30, 31] . It has also been used for studying the 
kaon interferometry in these collisions [32|. 



III. ANISOTROPIC FLOWS 

The anisotropic flows v n of particles are the Fourier 
coefficients in the decomposition of their transverse mo- 
mentum spectra in the azimuthal angle <fi with respect to 
the reaction plane |3^], i.e., 
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ll + J2 2v n{PT:y)cos{n<f>)} (2) 



Because of the symmetry cj> <-> — <f> in the collision ge- 
ometry, no sine terms appear in the above expansion. 
For particles at midrapidity in collisions with equal mass 
nuclei, which will be considered in the present study, 
anisotropic flows of odd orders vanish as a result of the 
symmetry (f> «-> 4> + n. The anisotropic flows generally de- 
pend on the particle transverse momentum and rapidity, 
and for a given rapidity the anisotropic flows at trans- 
verse momentum px can be evaluated from 



Vnipr) = (cos(n0)) 



(3) 



where (• • •) denotes average over the azimuthal distri- 
bution of particles with transverse momentum px- The 
anisotropic flows v n can further be expressed in terms of 
following single-particle averages: 
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where p x and p v are, respectively, projections of particle 
momentum in and perpendicular to the reaction plane. 

Since the AMPT model also provides information on 
the spatial anisotropy of colliding matter, which is re- 
sponsible for generating the momentum anisotropic flows, 
it is of interest to introduce the spatial anisotropic coeffi- 
cient s n by expressions similar to those for the anisotropic 
flows v n but in terms of the spatial distributions of par- 
ticles in the transverse plane. 
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FIG. 1: (Color online) Time evolutions of spatial anisotropic 
coefficients s n and anisotropic flows v n of partons in midra- 
pidity from Au + Au collisions at y/s = 200 AGeV and 6 = 8 
fm for parton scattering cross sections a p = 3 mb (left panels) 
and Up = 10 mb (right panels). 

Using the AMPT model in the string melting scenario 
with parton scattering cross sections a p — 3 and 10 mb, 
we have studied the time evolutions of s n and v n of par- 
tons in midrapidity from Au + Au collisions at y/s = 200 
AGeV and 6 = 8 fm, and they are shown in Fig. ^ 
The spatial anisotropic coefficients s n at a given time are 
evaluated using parton positions at their last interactions 
before this time as in Ref. 0] • I n this way, final values of 
s n are given by the parton spatial distributions at freeze- 
out. Similarly, final values of v n are determined by the 
parton momentum distributions at freeze-out, when they 
are converted to hadrons by coalescence. It is seen from 
Figs, ^a) and (b) that the spatial elliptic deformation 
S2 is initially large and decreases with time. It reaches 
the saturation value at about 10-15 fm/c for both par- 
ton scattering cross sections. The same behavior is ob- 
served for the fourth-order spatial anisotropic coefficient 
S4 except that it reaches the saturation value earlier at 
about 5 fm/c. The higher-order spatial anisotropic co- 
efficients sq and sg are essentially zero. The nonzero 
spatial anisotropic coefficients s n imply that the parton 



3 



spatial distribution is nonspherical at freeze-out. We note 
that the larger parton scattering cross section leads to a 
smaller spatial anisotropy at freeze-out. 

The time evolutions of the momentum anisotropic 
flows v n of partons are shown in Figs. ^Jc) an d (d) for 
the two parton scattering cross sections of 3 and 10 mb, 
respectively. The elliptic flow Vi is seen to saturate ear- 
lier in the collisions, i.e., at about 7 fm/c. This is similar 
for both parton scattering cross sections, except that the 
larger one leads to a larger elliptic flow. As to i>4 of 
partons, its value reaches the final value at a later time 
than that for v%, i.e., at about 10 fm/c, which may be 
due to the fact that the parton also affects its V4. 
For other higher-order anisotropic flows, Vq of partons is 
very small, having a magnitude less than 0.1% for both 
parton scattering cross sections, while vg of partons is 
essentially zero. An interesting result predicted by the 
AMPT model is that the parton V4 is more sensitive to 
the parton cross section than V2- For instance, the final 
values of v-i and V4 of partons are about 5.4% and 0.29% 
, respectively, for a p = 3 mb, and they are increased to 
about 7.5% and 0.67% , respectively, for er p = 10 mb. 
Therefore, V4 of partons is enhanced by a factor of about 
2.3 while their vi is only enhanced by a factor of about 
1.4, when the parton cross section is increased by a factor 
of about 3. 

The anisotropic flows of partons are transferred to 
those of hadrons when the latter are formed from the 
coalescence of quarks and/or antiquarks. Although scat- 
terings among hadrons are included in the AMPT model, 
they do not affect much the hadron anisotropy flows as 
a result of the small spatial anisotropy and pressure in 
the hadronic matter - In Fig. we show the fi- 
nal anisotropic flows v n of charged hadrons in the pseu- 
dorapidity range < 1.2 in minimum bias Au + Au 
collisions at s/s = 200 AGeV as functions of transverse 
momentum px for parton scattering cross sections a v = 3 
and 10 mb. Also shown in Fig. [21 are recent experimen- 
tal data from STAR collaboration [l4|. It is seen that 
the parton scattering cross section a p — 3 mb underesti- 
mates the data at higher p T (> 1 GeV/c) while a p = 10 
mb seems to give a better fit to the data. The values 
of hadronic vq are in agreement with the data within er- 
ror bars, although they are essentially zero. Again, U4 
of charged hadrons exhibits a stronger sensitivity to the 
parton cross section than their V2, especially at higher px- 
Therefore, V4 of hadrons provides a more sensitive probe 
of the initial partonic dynamics in relativistic heavy-ion 
collisions than their V2- 

We note that the parton scattering cross section used 
in the present study to reproduce the measured hadron 
elliptic flow is much smaller than that needed in Ref. [34[ 
based only on the parton cascade model. This is due to 
the fact that in Ref.[34| the measured hadron rapidity 
distribution and flow were compared to those of partons. 
Since the number of partons in our model is about twice 
that in Ref.[34[! we obtain a larger effect on parton el- 
liptic flow from parton scatterings with the same cross 
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FIG. 2: (Color online) Anisotropic flows V2 (a), V4 (b), and ve 
(c) of charged hadrons in the pseudorapidity range \rj\ < 1.2 
from minimum bias Au + Au collisions at y/s = 200 AGeV as 
functions of transverse momentum pr for parton scattering 
cross sections a p — 3 (open squares) and 10 (solid squares) 
mb. The experimental data are from STAR Collaboration 
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section. Furthermore, the elliptic flow of hadrons that 
are produced from parton coalescence is larger than that 
of partons if the increase of the parton elliptic flow with 
transverse momentum is less than linear, which is indeed 
the case as shown below in Fig[31 



IV. SCALINGS IN ANISOTROPIC FLOWS 

The experimental data indicate that there is a scaling 
relation among hadron anisotropic flows, i.e., v n (px) ~ 

v 2^ 2 (Pt) [13 • It has been shown by Kolb 35] that such 
scaling relation follows naturally from a naive quark coa- 
lescence model |36( that only allows quarks with equal 
momentum to form a hadron. Denoting the meson 
anisotropic flows by v nt M(pr), he found that i>4,m(pt) = 
(l/4)u| m(.Pt) h" quarks have no higher-order anisotropic 
flows. Since mesons dominate the yield of charged par- 
ticles in experimental data, the smaller scaling factor of 
1/4 than the empirical value of about 1 indicates that 
higher-order quark anisotropic flows cannot be neglected. 
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Including the latter contribution, one can show that 
w 4: m(2pt) ^ 1 lv 4i g(p T ) 



v IaA 2 Pt) 

V6,m(2pt) 
(2pr) 



J 2,M 



4 2«I i ,( Pr )' 

1 ( VA, q {PT) V6,q{pT) 

4 Ul^Cpr) «2, g (Pr) 



(8) 
(9) 



where v n ^ q {pT) denotes the quark anisotropic flows. The 
meson anisotropic flows thus satisfy the scaling relations 
if the quark anisotropic flows also satisfy such relations. 



The resulting hadron scaling factors of 3/4 and 1/2 are, 
however, smaller than the one extracted from measured 
anisotropic flows of charged hadrons. Since the naive 
quark coalescence model does not allow hadron formation 
from quarks with different momenta as in more realistic 
quark coalescence models p^. 12^.123 ? ], it is not expected 
to give a quantitative description of the experimental ob- 
servation. Such effects are, nevertheless, included in the 
AMPT model, which have been shown in Fig. [5] to re- 
produce the measured hadron anisotropic flows. 
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FIG. 3: (Color online) Transverse momentum dependence 
of midrapidity parton anisotropic flows V2, «4 and v§ from 
minimum bias events for Au + Au at i/s — 200 AGeV with 
parton scattering cross section 10 mb. Also plotted are v\ 
(open triangles) and v\ (open stars). 

In Fig. |21 we show the pt dependence of the 
anisotropic flows «2, t>4 and vq of midrapidity partons 
from minimum bias events in Au + Au collisions at 
\fs = 200 AGeV for a parton scattering cross section of 
10 mb. Also shown in Fig. [3] are v\ (open triangles) and 
v\ (open stars). It is seen that the parton anisotropic 
flows from the AMPT model indeed satisfy the scaling 
relation v nt q(pr) ~ ^ g (Pr)- With a parton scaling fac- 
tor of about 1, we obtain the following scaling relations 
among meson anisotropic flows: 

V4,m(pt) ^ 3 v 6 ,m{pt) ^ 1 
4m(?t) ~ 4' vl M {p T ) ~ 2' 



V. SUMMARY 

In summary, using the AMPT model in the string 
melting scenario, we have studied the anisotropic flows 
«2, i>4, and vq of both partons and hadrons in heavy- 
ion collisions at RHIC. We find that the anisotropic 
flows are developed early in the collision when the mat- 
ter is in the partonic stage and has appreciable spatial 
anisotropy. With a parton scattering cross section of 
about 10 mb, recent data on hadron isotropic flows from 
Au + Au collisions at y/s — 200 AGeV can be approxi- 
mately reproduced. Our results further show that i>4 is 
a more sensitive probe of the initial partonic dynamics 
in these collisions than t>2, and that higher-order par- 
ton anisotropic flows are non-negligible and satisfy the 
scaling relation v n , q (pT) ~ v 2qipr)- Using the quark 
coalescence model, we then show that the parton scaling 
relation leads naturally to the observed similar scaling re- 
lation among hadron anisotropic flows. Studies of higher- 
order anisotropic flows in relativistic heavy-ion collisions 
are thus expected to provide useful information on the 
collision dynamics and the properties of the hot dense 
matter during the initial stage. 
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